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Using a variational method we study a sequence of the one-electron atomic and molecular- 
type systems H, , H^ + ^ and //| 3+ ' in the presence of a homogeneous magnetic field ranging 
B = - 4.414 x 10 13 G. These systems are taken as a linear configuration aligned with the magnetic 
lines. For -ff^ 2+ ' the potential energy surface has a minimum for B ~ 10 11 G which deepens with 
growth of the magnetic field strength (JETP Lett. 69, 844 (1999)); for B > 10 12 G the minimum 
of the potential energy surface becomes sufficiently deep to have longitudinal vibrational state. We 
demonstrate that for the (ppppe) system the potential energy surface at B > 4.414 x 10 13 G develops 
a minimum, indicating the possible existence of exotic molecular ion Hi T \ We find that for almost 
all accessible magnetic fields is the most bound one-electron linear system while for magnetic 
fields B > 10 13 G the molecular ion H., becomes the most bound. 

I. INTRODUCTION 

In the presence of a strong magnetic field the properties of atomic and molecular systems change substantially 
leading to new unexpected phenomena. It has been a long time since the predictions made by Kadomtsev and 
Kudryavtsev jjj , and Ruderman , about the possible existence of unusual chemical compounds in the presence of a 
strong magnetic field which are aligned along the magnetic line. This phenomenon can be explained by the effective 
quasi-one-dimensionality of Coulombic systems appearing due to the enormous Lorentz force acting on a system in 
transverse direction. The transverse size of such systems is governed by the magnetic field and is of the order of the 
cyclotron radius p ~ ^~ x / 2 while the longitudinal size remains of atomic (molecular) order shrinking logarithmically 
~ l/log7 (see, for example, |J)Q. As a consequence of such quasi-one-dimensionality the electron- nuclei attraction 
becomes more effective, yielding the possibility of compensating the Coulombic repulsion of nuclei and thus to the 
possible existence of exotic systems that do not exist without a magnetic field. In particular, it was conjectured that 
the presence of a strong magnetic field can lead to the formation of linear molecules (linear chains) situated along 
magnetic lines [^|,^[j6) . Recently, the first quantitative study of the possible existence of the exotic molecular ion H^ 2+S> 
or, in other words, a bound state in the system of three protons and an electron (pppe) was carried out 0. It turned 
out that for B > 10 11 G the potential energy surface as a function of internuclear distances possesses a minimum, 
which deepens with magnetic field growth. It provides an evidence that for magnetic fields B > 10 11 G such molecular 
ion can exist 0. In the present article we extend this study and show that for magnetic fields B > 10 13 G a bound 
state in {ppppe) system or, in other words, another exotic molecular ion H\ can exist. 

The aim of this article is to perform a detailed quantitative study and comparison of the one-electron systems H, 
H^, H^ 2+ ^ and h[ 3+ ^ in strong magnetic field. Our study is limited to an exploration of the ground state only. It 
is assumed that the Born-Oppenheimer approximation holds, which implies that the positions of protons are kept 
fixed. We consider the linear molecular configuration for which the protons are aligned with the magnetic field (linear 
configuration parallel to a magnetic line). This configuration looks optimal from a physical point of view. Spin effects 
(linear Zeeman effect) are neglected. The magnetic field ranges from up to 4.414 x 10 13 G, where non-relativistic 
considerations hold (for a detailed discussion see, for instance, || and references therein). In order to study the 

(2+) 

stability of the system , the longitudinal vibrational modes of the nuclei arc considered, while transversal modes 

(deviation from linearity of the configuration as well as alignment along magnetic line) are not taken into account. 
The position of the lowest longitudinal vibrational state is estimated. 

Our calculations demonstrate that for all studied one-electron systems the binding energy = 7 — E to tai (affinity 
to keep the electron bound) grows with the increase of the magnetic field strength for B = — 4.414 x 10 13 G. For 
almost all accessible magnetic fields B < 10 13 G the molecular ion H% (see the results in jg]) is the most bound 

one-electron system while for magnetic fields B > 10 13 G the molecular ion H^ 2+ ^ becomes the most bound. Our 



1 Here 7 = B/Bo, where Bo = 2.3505 x 10 9 G, is in Ry and, in particular, it has a meaning of the energy of free electron in 
magnetic field. 
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calculations show as well that for magnetic fields B > 10 G the potential energy surface of the system (ppppe) 

(3+) 

exhibits the appearance of a minimum indicating the possible formation of such a molecular ion H y A '. It hints at 

(2+) 

a certain hierarchy of the thresholds for the appearance of new one-electron linear molecular system: iJg - at 
B > 10 n G, H[ 3+) - at B > 10 13 G etc Q. If this is so, it is very unlikely that the system H^ +) can appear for 
magnetic fields existing on the surface of neutron stars. For larger magnetic fields a reliable study requires taking 
into account relativistic corrections, which is beyond of the scope of this article. So, the question about the possible 
existence of H^ 4+ ^ is open. 

II. METHOD 

The present calculations are carried out in the framework of a variational method using, for each system, a unique 
simple trial function equally applicable to any value of the magnetic field strength. For a successful study a wise 
choice of trial functions is needed. A constructive criterion for an adequate choice of trial function (see (|] — (l^]) is 
the following: (i) the trial function ^t(x) should include all symmetry properties of the problem in hand; (ii) if the 
ground state is studied, the trial function should not vanish inside the domain where the problem is defined; (iii) the 
potential Vt (x) = , for which the trial function is an exact eigenfunction, should reproduce the original potential 
behavior near singularities as well as its asymptotic behavior. 

Let us first introduce notation (see Fig. p. We consider identical, infinitely- heavy centers of unit charge situated 
on the z-axis. The magnetic field of strength B is directed along the z axis, B = (0, 0, B). Throughout the paper the 
Rydberg is used as the energy unit. For the other quantities standard atomic units are used. 

The potentials corresponding to the systems we study are given by 
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where the quantity p = \J x 2 + y 2 is the distance from the electron to the z-axis, and the fj are the distances from 
the electron to the zth center. 

III. TRIAL WAVE FUNCTIONS 

A class of trial functions which satisfy the above-mentioned criterion incorporating the basic physical features of 
Coulomb systems in a magnetic field directed along the z-axis contains a function 

IpO = e -Vfri.".rN,p) ) 
TV 

<p(n,...r N ,p) =^a l r. i +/?5p 2 /4 (2) 
»=l 

as a basic element, where and /3 are variational parameters, and N is the number of centers. This function is nothing 
but the product of ground-state Coulomb functions and the lowest Landau orbital and is the exact eigenfunction of 
the potential 



2 In framework of non-relativistic consideration it can be easily checked that there exist magnetic fields for which a potential 
energy surface for h!^ + ^ etc. can develop a minimum. 
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J „2 Bo 

-2^a l --f3B + ^a l a ] {h l -h J ) + ^a l f3Bf-+l3 2 — P — , (3) 



hi ■ hj = — !— [p 2 + (z - z l )(z - Zj)] , 

r i r j 

where hi is unit vector in ith direction. This potential reproduces Coulomb-like behavior near the centers and 
two-dimensional oscillator behavior in the (x,y) plane at large distances (cf. (1) ). 

Since the centers are identical the problem possesses permutation symmetry. Hence the ground state trial function 
must be permutationally-symmetric with respect to permutation of positions of the centers, <-> rj. 

One of the simplest functions satisfying the above-mentioned criterion is the function (||) itself with the same a's 
(those can be placed equal to ai): 

^^-^(Et/fftV^ (4) 

where a%,^\ are variational parameters. Actually, this function is a Heitler-London type function multiplied by the 
lowest Landau orbital. If the intcrnuclcar distances are taken as parameters, it has in total (N + 1) variational 
parameters. It is quite natural from a physical viewpoint to assume that this function gives a description of the 
system near the equilibrium position. The potential V\{x) = , corresponding to this function is: 



2 n2 -i 



V, = Nal-Bp 1 -2a 1 J2[-)+^r- + ^lJ2^- fl i+ £ ( - ) . (5) 





(cf- (3)). 

Another possible trial function is a permutationally-symmetric linear superposition of the functions (0) in which 
all a's are placed to be equal to zero except one, which is denoted a-i 



N 

_ (6) 



i=l 



Here a2,/?2 are variational parameters. This function is a Hund-Mulliken type function multiplied by the lowest 
Landau orbital. One can assume that when all internuclear distances are sufficiently large this function dominates, 
giving an essential contribution. Thus, it describes an interaction of the hydrogen atom with (N — 1) charged centers. 
This function can also describe a possible decay mode of the studied system to the hydrogen atom and (N — 1) 
protons. As in Eq.(||), the trial function (||) is characterized by (N + 1) variational parameters. 

Another trial function describes a possible decay mode H^ N — > H% + (N — 2)p, N > 3, and is a 

permutationally-symmetric linear superposition of the functions (^|) with all a's are equal to zero except two, which 
are taken to be equal and denoted as a^: 

*3 = ( J2 e- a ^ +r ^e-^ B P 2 ' A , (7) 

i<j 

where 03,^3 are variational parameters. It seems that the function (^) gives the dominant contribution to the large 
internuclear distances when the distance between two neighboring centers is kept fixed. It can be true for magnetic 
fields B < 10 13 G, where the system with the lowest total energy is . Eq.(Q) also depends on (N + 1) variational 
parameters. 

Finally, for magnetic fields B > 10 11 G where the exotic state H^ 2 begins to appear we also introduce a trial 
function of the form 

*4=( £ e " Q4(r ' +rj+ri) je^ 4B ' ,2/4 . (8) 

i<j<k 

where again and (3^ are variational parameters. This function depends on N + 1 variational parameters. 

In order to incorporate in a single trial function the behavior of the system both near equilibrium and at large 
distances, we use an appropriate interpolation of the trial functions (Q), (|^), (f7j) and (|8|). There are several natural 
types of such interpolations: 
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(i) Total non-linear superposition of the form: 



5—nls — 



-a, r%—a 2 r 2 - 



{a' 1 ,...a' } 



(9) 



where a[, . . . a' N and (3' are variational parameters, and the sum is over all permutations of the parameters 
{a[, . . . a' N }. The function (^) is an A^-center modification of the Guillemin-Zener type function used for 
the description of the molecular ion H 2 in a magnetic field (see Eq.(2.5) in Q). If all parameters coincide 
a'i = ai, i = 1...N, the function (||) reduces to the Heitler-London type function (JI]). If the only parameter 
is non-zero, say a'y — ct2, this function reduces to the Hund-Mulliken type wave function (^), if two parameters 
are non-zero, and equal, say, oi-y = a' 2 = a% it reduces to the trial function ([?]), and if three parameters are 
non-zero, oi-y — a 2 = a' 3 = it reduces to (||). In general, there are 2N variational parameters (including the 
internuclear distances as parameters) characterizing the trial function (||) . 

(ii) Partial non-linear superpositions: 

these are special cases of (^|) when two or more parameters are equal. For instance, for N ~ 3 a\ — a' 2 ^ a' 3 : 

\ff 6 _ nls _ p = ^ e - a 'i( r i+ r 2)-a' 3 r 3 _|_ e -a' 1 (r 1 +r 3 )-a' 3 r 2 _|_ g-a'j (r 2 +r 3 )-at' 3 ri ~j p -P'Bp 2 /4 ^ 

This function can be considered as a non-linear interpolation between Eqs. (||) and (Q) for the H^ 2+ ^ ion. 

(iii) Linear superposition of Eqs. (Q), (||), (0), @ 

^ 7 -l s = Ax$i +A 2 * 2 +^3*3 +^4*4 , (11) 

where A\ 1 A2, A3, Ai are taken as extra variational parameters. It should be mentioned that in general ^1,2,3,4 
are not orthogonal and thus the parameters ^1^2,3,4 do not have a meaning of the weight factors. 

So, each of the above interpolations has a certain physical meaning describing different physical situations. In 
order to make up a rich description of the systems one can take a linear superposition of all previously mentioned 
functions with factors in front of them, those are treated as extra variational parameters. In this work we use only 
some particular cases of this general interpolation procedure (see below for details). 

The minimization procedure is carried out using the standard minimization package MINUIT from CERN-LIB. We 
use the integration routine D01FCF from NAG-LIB. All integrals are calculated with relative accuracy > 10 -8 . 



IV. RESULTS 



The results are summarized in Table I. In general, as a common feature the energy dependences for the hydrogen 
atom H and molecular ions H% and H£ +) are characterized by the growth of the binding energy, as well as a decrease 

of the equilibrium distances (R for H 2 , and R\ + R2 for H^ + ^) with increase of magnetic field strengths []. It seems 
reasonable to conjecture that this behavior persists for magnetic fields B > 4.414 x 10 13 G for which a non-relativistic 
treatment is no longer applicable. In Table I L is the "natural" size of the system which is determined by the 
internuclear distances: L = R for H 2 , L = Ri + R 2 for H^ 2+ ^ and L = Ri + i? 4 for H^ + \ see Fig. 1. For H^ 2+ \ 
the energy of the lowest vibrational state is defined as 5E — ^fno + + -^frus-, where lo± are determined by curvatures 
near minimum of the potential well along the principal axes: 1Z+ = R\ + i?2, 1Z- — R2 — Ri, see Fig. 1. One should 
emphasize that with a magnetic field growth all studied systems become more and more bounded. Also for both H 2 

(2+) (3+) 

and H^> the energy of the lowest vibrational state increases. The vibrational state for H\ was not studied. 



3 It is known that in the limit of large magnetic fields the binding energy of electron in H, Hi H n systems as well as in 
heavy atoms behaves oc (log7) 2 and the equilibrium distances oc (log7) _1 (see, for instance, [[l],^,^|,|l6| and discussion therein). 
Similar estimates should hold for other one-electron linear systems. 
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A. //-Atom 



Recently, an extensive numerical analysis of different features of the hydrogen atom in constant magnetic field was 
presented in the book by Ruder et al . |13| . In our study for the ground state of the iJ-atom we use the trial function 
of the form (§) at N = 1 (see @— (lj) 

where a and /3 are variational parameters. A comparison with the results obtained by Lai et al. pj shows that this 
simple trial function gives a relative accuracy 1CP 4 in energy for magnetic fields B = — 4.414 x 10 G. The binding 
energy for this system increases with magnetic field growth and the longitudinal size reduces drastically from 1.5 a.u. 
at B = to 0.242 a.u. at B — 4.414 x 10 13 G. The total energy for the iJ-atom is larger than the total energy for 
H2 for all studied magnetic fields, as well as larger than the total energies for H^ 2+S> and for H^ 3 ^ for B > 10 11 G 
and B > 4.414 x 10 13 G , respectively (see Table I). So, the i/-atom is a least bound one-electron system among H^, 
andfff +) . 

B. H2 molecular ion 

For this system a linear superposition of trial functions (^), (||) and (^|) for N = 2 is used for present variational 
calculation (cf. Q). In total there are ten variational parameters and the present calculation provides the most 
accurate ground state total energy for B > 10 9 G. The ground state corresponds to the positive parity state 1<t u . 
For magnetic fields in the region B = — 10 13 G this system is characterized by the lowest total energy among the 
one-electron systems we consider, being the most stable. The binding energy increases (and the equilibrium distance 
R decreases) with the increase of the magnetic field. It is worth noting that the longitudinal localization length of 
the electron (|z|) is larger than the "natural" size of the system L = R for magnetic fields B > (See Figs. 2,3). 
The energy of the lowest longitudinal vibrational level is estimated using a harmonic approximation of the bottom of 
the potential energy curve. The position of the lowest vibrational state in comparison with the Born-Oppcnhcimcr 
minimum of the total energy grows steadily from 5E = 0.01 Ry at B = to SE = 1.23 Ry for B = 4.414 x 10 13 G, 
respectively. 



C. H$ molecular ion 

For this exotic system H^ 2+ ^ theoretical evidence was provided of its existence for magnetic fields B > 10 11 G [Q. 
For the present variational calculations we use a linear superposition of trial functions ([|), (|[), (Q) and (Q) for N = 3. 
We have in total fourteen variational parameters. As in the case of i?^ the binding energy increases, and the natural 
size of the system L = R\ + R2 decreases, as the magnetic field grows. In contrast to the ion, the longitudinal 
localization length of the electron ( | z | ) is smaller than the natural size of the system. In this case the electronic cloud 
mainly surrounds the central proton (See Figs. 2,3). The potential energy surface on the plane (Rx^R-i) exhibits a 
well-pronounced minimum at R\ = R2 starting from B ~ 10 11 G. The energy of the lowest longitudinal vibrational 
state is calculated using a harmonic approximation around the minimum along the principal axes 1Z+ = R\ + R2, 
1Z- = i?2 — Ri in the plane (i?i,i?2)- Our results show that for B — 10 11 G the energy of the lowest longitudinal 
vibrational state in harmonic approximation SE = 0.189 Ry lies much above the top of the barrier AE = 0.054 Ry. 
However, the situation changes with growth of magnetic field. For B > 10 12 G the potential well becomes sufficiently 

deep to keep the lowest vibrational state inside of the well (see Table I). For B = 10 11 — 10 13 G the system H^ 2+ ^ is 

unstable towards a decay to +p. With the growth of magnetic field the total energy of H^ 2+ ^ becomes lower than 
the total energy for ■ This happens at B « 4.414 x 10 13 G and it means that for B > 4.414 x 10 13 G the exotic 

(2+) 

molecular ion becomes the most bound one-electron system. 



D. H\ molecular ion 

To study this system we use a superposition of the trial functions (Q), (Q), (Q) and (||) for N = 4 in total with 
fifteen variational parameters. A symmetry argument implies that the minimum, if it exists, should correspond to the 
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configuration R3 = R2 and R4 = R\ (see Fig. |]). The results show that for magnetic fields B < 10 13 G the potential 
energy surface in the plane {R\, R2) does not exhibit any minimum. The first indication to a possible existence of the 
minimum comes at B = 10 13 G and for B = 4.414 x 10 13 G the potential energy surface indicates the appearance of a 
well-pronounced minimum of depth AE — QARy at R\ — 0.225 and R2 — 0.058. The longitudinal localization length 
of the electron (\z\) is smaller than the natural size of the system L = 2Ri similarly to what takes place for H% . 
The electronic cloud surrounds mainly the two central protons (see Figs. 2,3) similar to ' . It is important to 
stress that the total energy for this system is smaller than the total energy of the H atom but is larger than the total 
energy for H+ and Hf +) at B = 4.414 x 10 13 G. 



V. CONCLUSIONS 



In the presence of a strong magnetic field the formation of exotic chemical compounds can occur. In particular, 
in present article we demonstrate that in addition to the ii-atom and the B^ molecular ion two more one-electron 

(2+) (3+) 

systems B^ and H\ can exist which do not exist without a magnetic field. The results of our consideration show 
an appearance of a minimum in the potential energy surfaces for the systems (pppe) and (ppppe) for magnetic fields 
B > 10 11 G and at B = 4.414 x 10 13 G, respectively, indicating the formation of the exotic ions H^ 2+ ^ and H^ + \ 
Among these systems the molecular ion possesses the lowest total energy in the region B < 10 13 G being the 

most stable towards a fission of proton(s) or an electron. The ion B^ 2+ ^ can decay to H% + p but not to B + p + p. 
However, for magnetic fields B > 4.414 x 10 13 G the ion B^^ becomes the most stable one-electron system possessing 

(24-) 

the lowest total energy. Our estimates show that the lowest longitudinal vibrational state of B^ can also exist for 
magnetic fields B > 10 12 G. 

For a magnetic field B ~ 4.414 x 10 13 G, where we assume that a non-relativistic treatment is still valid, the exotic 
molecular ion can exist. This ion is unstable towards a decay to B^ +p + p and B^ +S> +p, however having no 

decay mode B + p + p + p. A further study of this system beyond the non-relativistic treatment is needed for reliable 
evidence of its existence. We should also mention that there is a general question about applicability of the Born- 
Oppenheimer approximation for the systems in a magnetic field (see discussion in f p0] , [i"3"p |] ) . In all our calculations 
we use the standard Born-Oppenheimer approximation. 
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FIG. 1. One-electron systems in a magnetic field B - explanation of the notations used, (a) H-aXom, (b) , (c) , (d) 

i?4 . The positions of charged centers are marked by bullets. 
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FIG. 2. Electronic probability density for (a) # 2 + , (b) h£ +) and (c) Hf +) for B = 4.414 x 10 13 G. Normalization of |*| 2 is 
not fixed. 
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FIG. 3. Electronic probability density contours for (a) H% , (b) Hf +) and (c) Hf +) for B = 4.414 x 10 13 G . The position 
of the centers is indicated by small circles. The longitudinal localization length of the electron (\z\) is displayed. 
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TABLE I. Variational results for the one-electron linear systems H, H£ , and H 4 ina strong magnetic field for total 

energy E, the natural size of the system L (see text), longitudinal electron localization length (\z\), longitudinal vibrational 
energy SE measured from the bottom of the well, and the height of the barrier AE (see text for further explanations). Results 
in parenthesis are taken from * Lai et al. ||, ** Lai et al. ||, * Wind * Wille Q, ** Kravchenko and Liberman ° 

Z I ■ 



B (Gauss) 




H 


Ht 








E 


-1.000 


-1.20525° 










— 


(-1.205268+) 


— 


— 




L 


— 


1.997 


— 


— 






— 


(2.000*) 


— 


— 





(N> 


1.50 


1.704 


— 


— 




SE 


— 


0.0107 


— 


— 






— 


(0.0106 n ) 


— 


— 




AE 


— 


— 


— 


— 




E 


36.996 


35.036° 


36.429°° 


— 






(36.929*) 


(35.043*) 


— 


— 




L 


— 


0.593 


1.606 


— 






— 


(0.593*) 


— 


— 


10 11 


(1*1) 


0.685 


0.624 


0.864 


— 




SE 


— 


0.084 


0.189 


— 






— 


(0.074**) 


— 


— 




AE 


— 


— 


0.054 


— 




E 


413.94 


408.30° 


410.296°° 


— 






(413.57*) 


(408.56*) 


— 


— 




L 


— 


0.283 


0.692 


— 






— 


(0.278*) 





— 


10 12 


(N) 


0.441 


0.348 


0.438 


— 




SE 


— 


0.259 


1.033 


— 






— 


(0.243 **) 


— 


— 




AE 


— 


— 


0.833 


— 




E 


4232.6 


4218.665° 


4220.090°° 








(4231.6*) 


(4218.674*) 


— 


— 




L 





0.147 


0.330 











(0.15*) 






10 13 


(1*1) 


0.301 


0.214 


0.242 






SE 




0.714 


2.406 










(0.625**) 








AE 






3.58 






E 


18747.4 


18724.501 


18723.88 


18735.70 




L 




0.102 


0.220 


0.450 


4.414 x 10 13 


(1*1) 


0.242 


0.164 


0.168 
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SE 




1.233 
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